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Abstract 

Dark Matter (DM), arising from an Inert Higgs Doublet, may either be light, below the 
W mass, or heavy, above about 525 GeV. While the light region may soon be excluded, 
the heavy region is known to be very difficult to probe with either Direct Detection 
(DD) experiments or the Large Hadron Collider (LHC). We show that adding a second 
Inert Higgs Doublet helps to make the heavy DM region accessible to both DD and 
the LHC, by either increasing its couplings to the observed Higgs boson, or lowering its 
mass to 360 GeV < moM, or both. 
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1 Introduction 


The long-awaited Higgs boson, with a mass of rrih ~ 125 GeV, was famously discovered in 
2012 by the ATLAS and CMS experiments at the CERN Large Hadron Collider (LHC) [I112]- 
Although its properties are in accordance with the predictions of the Standard Model (SM), 
including Electro-Weak (EW) precision data, it remains an intruiging possibility that the 
observed Higgs boson, denoted here as h, may just be one member of an extended Higgs 
sector. A good motivation for the latter is the idea that it might provide a candidate for 
Cold Dark Matter (CDM). 

Although the nature of Dark Matter (DM) is not yet known, according to the Standard 
Cosmological Lambda-CDM Model [3] it should be a particle which is stable on cosmological 
time scales, cold, i.e., non-relativistic at the onset of galaxy formation, non-baryonic, neutral 
and weakly interacting. Various candidates for such a state exist in the literature, the most 
well-studied being the Weakly Interacting Massive Particles (WIMPs) [H El E] , with masses 
between a few GeV and a few TeV. Any such WIMP candidate must be cosmologically stable, 
usually due to the conservation of a discrete symmetry, and must freeze-out (i.e., drop out 
of thermal equilibrium) to yield the observed relic density [SQ 

= 0.1199 ±0.0027. (1) 

It is clear that the SM Higgs sector cannot provide a WIMP candidate, since its Higgs 
boson is unstable. However, it was suggested some time ago that the Higgs sector could be 
extended by the addition of an extra doublet, which may not develop a Vacuum Expectation 
Value (VEV), leaving a discrete Z 2 symmetry unbroken [7j. Independently, it was later 
shown that an extra scalar doublet with zero VEV, odd under a discrete Z 2 symmetry, could 
yield monojets at hadron colliders while being constrained by DM considerations (the hrst 
time to our knowledge that any connection with hadron colliders or DM was made) [8]. 
This possibility, which became known as the Inert Doublet Model (IDM), has been studied 
extensively for the last few years (see, e.g., p uni EH). Since the IDM involves 1 Inert 
Doublet plus 1 active Higgs Doublet, we shall also refer to it henceforth as the I(1±1)HDM. 

In the IDM, aka the I(1±1)HDM, one extra spin-zero SU{2)l doublet with the same 
quantum numbers as the SM Higgs doublet is introduced. One of the possible vacuum states 
in this model involves the first doublet acquiring a VEV, henceforth called the active doublet, 
while the second doublet does not develop a VEV and is referred to as the inert doublet since 
it does not take part in EW Symmetry Breaking (EWSB). Since this doublet does not couple 
to fermions and it is by construction the only Z 2 -odd held in the model, it provides a stable 
DM candidate, namely the lightest state among scalar and pseudo-scalar Z 2 -odd particles. 

The I(1±1)HDM remains a viable model for a scalar DM candidate, being in agreement 
with current experimental constraints. As of now, there are two regions of DM masses where 

^Since the Planck 2015 results quotes various results for depending on which spectra such as TT, 

TE and EE are used, we prefer to use here the Planck 2013 result whose error encompasses all of them. 
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one can expect viable solntions: a low DM mass region, 53 GeV < rriDM ^ and a heavy 
DM mass region, moM ^ 525 GeV. The most recent experimental data, both from direct 
detection experiments and from the LHG, has rednced the viable parameter space in the low 
mass region [HI OS]. In the heavy mass region, however, where the sensitivity of DM direct 
detection experiments decreases signihcantly with increasing DM mass, the DM candidate 
may escape possible detection in this model. 

In a recent paper na we studied DM in a model with 2 inert Higgs plus 1 active Higgs 
doublet, which we referred to as the I(2+1)HDM. In particular we focused on the region of 
parameter space of the I(2+1)HDM where the DM candidate, the lightest inert scalar, is in 
the light mass region {rriDM ^ mw)- We found that the extended scalar sector can relax 
the exclusion limits from direct detection experiments, providing a viable DM candidate in 
a region of parameter space which would be excluded in the I(1+1)HDM. In this paper we 
study the heavy DM mass region of the I(2+1)HDM. We show that heavy Higgs DM in 
this model becomes more readily observable as a result of either lowering the DM mass to 
360 GeV < ttidm, or increasing the DM-Higgs coupling, or both, while always maintaining 
the DM relic density within the required region. 

The layout of the remainder of this paper is as follows. In section we review the 
I(2+1)HDM and focus on a simplihed version of the model based on a smaller number of 
parameters. In section we calculate the relic density in the I(2+1)HDM, discussing the 
relevant DM annihilation scenarios, including the extended co-annhilating (pseudo-)scalar 
sector. Section will be focused on new features of the I(2+1)HDM with respect to the 
I(1+1)HDM in the context of DM phenomenology, including enhanced DM-Higgs couplings 
and the new mass region 360 GeV < ttidm ^ 525 GeV as well as on discussing the resulting 
improved prospects for direct detection. In section]^ we present heavy DM signals via Higgs 
mediation at the LHG in the I(2-|-1)HDM which look more promising than in the I(1-|-1)HDM. 
Finally in section]^ we conclude the paper. 

2 The I(2-|-1)HDM 

2.1 The scalar potential 

It has been shown in [T^ that an N-Higgs-Doublet Model potential symmetric under a group 
G of phase rotations can be divided into two parts; a phase invariant part, Vq, and a collection 
of extra terms ensuring the symmetry group G, Vg- 

We construct our Z 2 -symmetric 3-Higgs Doublet Model potential generated by the group 

^ = diag(-l,-l,l). (2) 
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which is of the following forn 0 

Vi(2+i)hdm = Vo + Vz2 (3) 

V^O = ~/^l(010l) ~ A''2(*?^202) ~/^3(0303) 

+ All(0j0l)^ + A22(0202)^ + ^33(0303)^ 

+ Ai2(0l0l)(0202) + A23(02‘/’2)(0303) + A31 (0303) 

+ Ai2(0l02)(020l) + A23(02</>3)(0302) + A31 (0^01) (0l03) • 

Vz2 = ~/^12(0102) + Ai( 0302)^ + A2(0203)^ + A3(030l)^ + h.C. 

We shall not consider CP-violation in this paper, therefore, we require all parameters of 
the potential to be real. 

The doublets are defined as 


01 — I H^+iA\ ) , 02 — I H°+iAl ) , 03 — [ v+h+iG° ) i (4) 

V J V ) \ V-2 J 

where 0 i and 02 are the two inert doublets and 03 is the one active doublet which plays the 

role of the SM Higgs doublet, with h being the SM-Higgs boson and G^, are the would-be 
Goldstone bosons. 

We assign Z 2 charges to each doublet according to the Z 2 generator in Eq.([^: odd-Z 2 
charge to the inert doublets, 0i and 02 , and even-Z 2 charge to the active doublet, 03. It 
is clear that the symmetry of the potential is respected by the vacuum alignment ( 0 , 0 , 

The neutral fields from the inert doublets could then in principle be DM candidates. These 
neutral fields are stabilised from decaying into SM particles as a result of the conserved Z 2 
symmetry of the potential after EWSB. 

To make sure that the entire Lagrangian and not only the scalar potential is Z 2 sym¬ 
metric, we assign an even Z 2 parity to all SM particles, identical to the Z 2 parity of the 
only doublet that couples to them, i.e., the active doublet 03. With this parity assignment 
Flavour Changing Neutral Currents (FCNCs) are avoided as the extra doublets are forbidden 
to couple to fermions by Z 2 conservation. 

The Yukawa Lagrangian of the model is identical to the SM Yukawa Lagrangian, with 03 
playing the role of the SM Higgs doublet: 

^Yukawa — ^ rnnQiri,L<p3Un,R + ^mnQrri, L4>3dn,R 

+^mrJm,L4>3en,R + ^mJm,L^3^n,R + h.C. (5) 

^Note that adding extra Z 2 -respecting terms such as (030i)(0j03), (0|02)(0303)) (0i02)(0i0i) and/or 
( 01 02)(0202) does not change the phenomenology of the model. The coefficients of these terms, therefore, 
have been set to zero for simplicity. 
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2.2 Mass eigenstates 

The minimum of the potential sits at the point (0, 0, with 


Mi 

-^33 ‘ 


The mass spectrum of the scalar particles are as follows. 

• The fields from the active doublet 

The fields from the third doublet, G^,G^,h, which play the role of the SM Higgs 
doublet helds have squared masses: 

uIqo = rnQ± = 0 

ml = 2iil (6) 


The CP-even neutral inert fields 

The pair of inert neutral scalar gauge eigenstates, H 2 , which are rotated by 


Re,, — 


cos 9h sin 9h 
— sin 9h cos 9h 


with tan26*/i = 


2/i 


12 


hi ~ R<i>i ~ h2 + ^ 


</>2 


into the mass eigenstates, Hi,H 2 , have squared masses: 

^Hi = (“hi + ^01) cos^ + {-lA + ^02) - 2hi2 sin 9h cos 9h 

'^H 2 = (“hi + A0i) sin^ 9h + (-h2 + ^02) cos^ 9h + sin 9h cos 9h 

where = “(-^si + -^31 + 2A3)u^, A(^2 = -(A23 + A23 + 2A2)u^ 


(7) 


The charged inert fields 

The pair of inert charged gauge eigenstates, 0^, 0^, which are rotated by 


Re^ — 


cos 9 c sin 9 c 
— sin 9r cos 9c 


with tan26*r = 


2/1' 


12 


h? - RA - hi + a; 


02 


into the mass eigenstates, have squared masses: 


m]^± = {-lA + A'^J cos^ 9c + (-hi + - 2/1^2 sin 9c cos 9, 

nAj± = (-/ii + A'^J sin^ 9c + (-hi + R' 4 , 2 ) cos^ + 2,^12 sin 9c cos 9, 


where = ^(^3l)n^ R^ = ^(^23)^^ 


( 8 ) 


The CP-odd neutral inert fields 

The pair of inert pseudo-scalar gauge eigenstates, A°, A2, which are rotated by 


Re,. — 


cos 9 a sin 9 a 

— sin 9a cos 9a ' ’ 


with tan20„ = 


2/i' 


12 


hi R'A 9-2 + A 


// 

(p2 


4 



into the mass eigenstates, Ai, ^ 2 , have squared masses: 

= (-hi + A'^i) COS^ + (-h2 + ^ 02 ) “ 2/1i 2 6'a COS 

"^A2 = (-h? + ) sin^ Qa + (-h2 + A'^J cos^ + 2/i^2 sin Qa cos 

where = -(A 31 + Ag^ - 2 A 3 )n^, = ^('^23 + A 23 - 2 A 2 )n^ (9) 

We will refer to (i^i, Ai, H^) as the helds from the hrst generation and to {H 2 , A 2 , H^) 
as the helds from the second generation. Each of the four neutral particles could, in prin¬ 
ciple, be the DM candidate, provided it is lighter than the other neutral states. In what 
follows, without loss of generality, we assume the CP-evenj^ neutral particle Hi from the hrst 
generation to be lighter than all other inert particles, that is: 

rriH^ ( 10 ) 

(Note that this choice is arbitrary: if the CP-even particle from the second generation, H 2 , 
where to be assumed lighter than the other inert states, then H 2 will play the role of the DM 
candidate.) 

Assuming the CP-even neutral inert particles are lighter than the CP-odd and charged 
inert particles puts the following constraints on the parameters: 

2A2,2A3<A'23,A'3i<0. (11) 

In our DM analysis, we consider cases where the mass alignment is changed, but where Hi is 
always the lightest inert state and hence is the DM particle. In the remainder of the paper 
the notations Hi and DM particle will be used interchangeably. 

2.3 Simplified couplings in the I(2+1)HDM 

Due to the large number of free parameters in the I(2-|-1)HDM which makes it impractical to 
analyse the model in the general case, we focus on a simplihed case where parameters related 
to the hrst inert doublet are k times the parameters related to the second doublet 

A 3 = kX2, A 31 = k\23, Ag]^ = A:A23, ( 12 ) 

resulting in 

= kA^,, a;^ = a:a;^, a; = a;A';, ( 13 ) 

without introducing any new symmetry to the potential. The motivation for this simplihed 
scenario is that in the k = 0 limit the model reduces to the well-known I(1-|-1)HDM. We 

^For the CP-even particle to be the DM candidate rather than the CP-odd particle, it is required that 
niHi < rnAi j which leads to A2, A3 < 0 . If instead Ai is assumed to be the DM candidate, A2, A3 > 0 . Hence, 
the results of our analysis are also applicable to the Ai DM case by changing the sign of A2 and A3. 
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assume no specific relation among the other parameters of the potential. It is important 
to note that the remaining quartic parameters do not influencethe discussed DM and LHC 
phenomenology of the model and thus their values have been fixed in agreement with the 
theoretical constraints discussed in the original paper [T3] and compliant with the results on 
unitarity obtained in |16j . 

The k = 1 case 

In this paper we focus on the k = 1 case in the heavy DM mass regiorj^ The mass spectrum 
in this case is simplihed to: 

= -fA + ^2 - + 2/^L, (14) 

m%± = -/i^ + 

+ K2 - k-A, 

The quartic couplings in the potential can be written in terms of the masses of the physical 
particles as: 


^23 = 

^2 = (15) 

2 

A23 = gHiHih —^("^^1 ~ 

where gu^Hih = ^23 + A23 + 2A2 is the Higgs-DM coupling. The Feynman rules for this model 
are presented in Appendix 


3 Calculating the relic density in the I(2+1)HDM 

The relic density of the WIMP (identified in our model as the lightest inert scalar Hi) is 
calculated with the assumption that the WIMP was in thermal equilibrium with the SM 
particles after inflation. Once the rate of 

DM DM o SM SM 

reactions becomes smaller than the Hubble expansion rate of the Universe, the WIMP freezes 
out, i.e., drops out of the thermal equilibrium. After freeze-out the co-moving WIMP density 
remains essentially constant with the current value estimated by the Planck experiment to 
be the one already given in Eq. ([^. 

^Other scenarios with k ^ 1 are studied in [T4] . 
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As mentioned, in the I(2+1)HDM one of the neutral inert (pseudo-)scalar particles play 
the role of the DM. The relic density of a (pseudo-) scalar DM candidate, S, after freeze-out 
is given by the solution to the Boltzmann equation: 


dns 

dt 


-3Hns - {(Teffv){nl - 


S = H,,H2,A,,A2, 


(16) 


where the thermally averaged effective (co)annihilation cross-section contains all relevant 
scattering processes of any SiSj pair into SM particles: 




eq eq' 

n/ n/ 


where 


n. 


eq 


~ exp(- 


n 


rrii - ms 
T 


). 


(17) 


(18) 


Therefore, only processes in which the mass splitting between a state Si and the lightest Z 2 - 
odd particle S {Hi in our case) are comparable to the thermal bath temperature T provide 
a sizeable contribution to this sum. 

In the I(2-|-1)HDM, the presence of additional inert particles has important consequences 
in the heavy mass regime. For lighter masses the most important channel for the annihilation 
of DM particles is the Higgs-mediated process 


HiHi ^ // 


see Fig. 12 r), as studied in HI. However, coannhilation with H 2 , Ai and A 2 may change 


the results significantly (see Fig. 12 d). 

For heavier masses the diagrams including one or two virtual gauge bosons, shown in Fig. 
13 also contribute to the total annihilation cross-section. Finally, co-annihilation plays an 
important role in scenarios with multiple particles which are close in mass. This scenario 
is realised in the I(2-|-1)HDM for the heavy DM mass region. Particles up to 20% heavier 
than the DM candidate may influence the DM relic density. Therefore, the co-annihilation 
diagrams should be included in calculating the effective annihilation cross-section. These dia¬ 
grams are presented in Figs. I^andf^- representing pure gauge channels and coannhilation 


channels involving the SM-like Higgs particle, respectively. 


3.1 Relevant co-annihilation scenarios 

In the I(2-|-1)HDM, the strength and importance of coannhilation processes depend on the 
mass splittings between the inert particles. We define 6a and Sc as the splitting between Hi 
and the pseudoscalar and charged state from the first generation, respectively, 

6A = mA^-mH^, 6c = mj^± -mHi- (19) 
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6a, c are related to the quartic couplings in the potential, which are constrained by the 
perturbativity (and unitarity) conditions, i.e., the Ads cannot be too large. As a result of 
this, all particles within one generation will have a similar mass. These masses, however, 
could have high values because of the (almost) unconstrained quadratic parameters /i| and 
h'12- 


- /^12 + ( 20 ) 

It is important to stress that, even if bounds on Aj were relaxed leading to larger values 
of 6a,c, there exist very stringent limits from relic density analysis. Coannihilation must 
occur at least between Hi,Ai and to achieve DM relic density in agreement with the 
current experimental measurements. This is a pattern followed by all general heavy scalar 
DM models. In the absence of these co-annihilation channels, the maximum relic density 
that can be achieved through HiHi —)■ SM SM (even when HiHj —)• SM SM is allowed) is 
of order 10“^ which is well below the observed value. 

The other important mass splitting. A, is defined as the mass difference between Hi and 
the other CP-even state H 2 (“splitting between doublets”): 

A = mH2 - (21) 

A is related to the quadratic parameter /if2 through 

/^?2 = = \ + ‘^rnnA) ■ ( 22 ) 

Note that fii 2 is not limited by any theoretical constraints - similar to /i| - and therefore 
A can in principle be very large. So, unless A is forced to be small by limits put on 
should also consider a case where the second doublet is decoupled from the first, leading to 
a scenario which was not listed in [TT]. Therefore, in the very heavy mass region one can 
consider: 

• Case G: with small ^A, Sc, A, where all inert particles are close in mass and co¬ 
annihilate with each other. 

• Case H: with small 6a, 6c and large A, where the second generation is effectively 
decoupled from the first generation and does not influence relic density calculations. In 
this case, the relevant diagrams are the (co)annihilation channels between fields from 
the lighter generation only. Hi, Ai, H^. 

Tab. [^summarises the two possible scenarios for relic density studies. 
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Small 6a, 6c 

Large 6a, 6c 

Small A 

Gase G is realised 

Goannhilation between Hi, H 2 is not efficient 
enongh and DM density is below experimen¬ 
tal bounds. 

Large A 

Gase H is realised 

There are no co-annihilatin channels open 
and gange annihlation rednces DM relic 
density ehectively below the experimental 
bonnds. 


Table 1: Valid regions of the parameter space schematically shown in terms of 6a, 6c and A 
in the heavy mass regime of I(2+1)HDM. 


3.2 The gauge limit 

To illnstrate the difference between cases G and H, let ns hrst consider the gange limit in 
both scenarios, which is the limit where all qnartic conplings Aj are set to zero. Therefore, 
all scalar self-conplings, inclnding the DM-Higgs conpling, are removed in this limit. As a 
resnlt 6a, c = 0, leading to degenerate Hi, Ai, states. Note that this limit is exclnded 
by resnlts of direct detection experiments, nevertheless, it is an interesting limit to stndy as 
it represents the main difference between cases G and H. In this limit Hi annihilates solely 
throngh the gange annihilation channels presented in Fig. 


14 


Non-zero Aj will lift this degeneracy and, at the same time, rednce the effective annihila¬ 
tion cross-section for a given mass. Therefore, the gange limit corresponds to the minimnm 
valne of m^i, for which it is possible to obtain a proper relic density for any valne of the 
Higgs-DM conpling. 

These resnlts are presented in Fig. for the two scenarios: case G where all particles 
have degenerate mass (in the 6a,c, A —)■ 0 A* —)• 0, 0 limit) and case H with large 

A where the second generation is deconpled from the hrst generation and all particles in the 
hrst generation are degenerate in mass (in the 6a,c —)■ 0 A* —j- 0 limit). 

It is clear that, for a given mass of mHi the destrnctive interference between an increased 
nnmber of coannhilation diagrams in case G leads to a rednced cross-section, i.e. larger DM 
relic density with respect to case H. The important conseqnence of this interference is that 
for case G it is possible to obtain proper relic density for smaller masses of DM candidate in 
comparison to case H. Note also that case H behaves like the I(1-|-1)IIDM (the Inert Donblet 
Model) in this limit. This similarity in behavionr will be repeated as we will show in the 
following sections. 
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Gauge limit for the l(2 + l)HDM 



mHj [GeV] 

Figure 1: The gauge limit in case G (blue curve) where all particles have degenerate mass 
(in the Sa,c, A —0) and case H (red curve) with large A where the second generation is 
decoupled from the first generation and all particles in the first generation are degenerate in 
mass (in the 5a,c 0 limit). 


3.3 The benchmark points 


With non-zero scalar couplings and mass splitting more diagrams contribute to the co¬ 
annihilation of DM - all diagrams shown in Figs. 14 and 15 contribute to the total an¬ 


nihilation cross-section. Here we present two benchmarks points, two sets of parameters, for 
which we have studied the DM relic density: 


• For case G with 6a = ^ GeV and A = 1 GeV 

Here all inert particles have similar masses and therefore can co-annihilate with each 
other. The degeneracy between charged and “pseudo-scalar” particles is allowed and 
doesn’t lead to any unacceptable results. The important degeneracy which must be 
avoided is Hi-Ai degeneracy leading to the scattering through the Z boson which is 
tightly constrained by direct detection experiments and puts a lower limit on 6a- 

• For case H with 6a = 6c = 1 GeV and A = 100 GeV 

Here the second generation of inert scalars is signihcantly heavier than the hrst one. 
Within each generation, however, particles are almost degenerate. 


Note that there are certain differences between cases G and H. In case H, the heavier 
generation of inert particles is decoupled from the hrst generation particles and does not 
inhuence the relic density calculations. The model in this case resembles the I(1-|-1)HDM. 
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Furthermore, in case G the Higgs-DM couplings which result in a relic density in agreement 
with experiment are larger in comparison to case H for the same DM mass. This difference 
is explicit in Fig. 





CaseG 
CaseH 
Planck+3 cr 
Planck-3 cr 


Figure 2: Relic density plots in case G (A = 1 GeV, blue line) and H (A = 100 GeV, red 
line) for = 550 GeV. The dashed horizontal lines show the 3a relic density limits from 
Planck in Eq. Q. 

This leads to the fact that for case G we can obtain viable relic density values for rriDM 
much smaller than in case H (or the I(1+1)HDM) in which the minimal value of muM 
resulting in DM relic density in agreement with Planck limits is rriHi ~ 525 — 535 GeV. In 
case G (with A = 1 GeV), however, the DM mass can be as low as ~ 375 GeV. This result 
is shown in Fig. which represents relic density plots for cases G (left) and H (right). Note 
that in case G the minimum rriHi which touches the lowest acceptable relic density limit (the 
green solid line) is 375 GeV (for a given A of 1 GeV), whereas in case H this minimum value 
is 525 GeV (the solid red line). 

Fig. I^is meant to represent the same benchmark points as in Fig. in the mH^-gHiHih 
plane. The bands correspond to proper relic density in agreement with Planck measurements 
in case G (for an exemplary A = 1 GeV, 5 = 1 GeV) in red and case H (for an exemplary 
A = 100 GeV, 5 = 1 GeV) in red. Note that, for the a given DM mass (and same 5), the 
Higgs-DM coupling in case G is much larger than in case H. 

4 DM in the I(2+1)HDM: direct detection 

4.1 Changes in A, 6 a , 6 c 

For the two viable scenarios G and H we allow the physical parameters to vary in the following 
regions: 

100keV<5^,5c<15GeV. (23) 
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Case G 


niH, = 775 GeV 


niH, = 725 GeV 


Case H 



- mH, = 775 GeV 

-mH, = 725 GeV 

- mH, = 675 GeV 

- mH, = 625 GeV 

. mn, = 575 GeV 

mn, = 525 GeV 
- mn, = 500 GeV 


Figure 3: Relic density plots for case G (left) with A = 1 GeV, 5 = 1 GeV and case H (right) 
with A = 100 GeV, 5 = 1 GeV. Note the solid red line in the case H which represents the 
minimum DM mass rriHi — ^25 GeV which just touching the lower relic density limit. The 
relic density plot for the same DM mass has been highlighted in case G (red solid line) which 
is well within the acceptable relic density limits. The dashed horizontal lines show the 3a 
relic density limits from Planck in Eq. (Q. 


9h, H, h 



mH,[GeV] 


— case G 

— case H 


Figure 4: Relic density bands in agreement with Planck measurements in case G (for an 
exemplary A = 1 GeV, 5 = 1 GeV) in red and case H (for an exemplary A = 100 GeV, 5 = 1 
GeV) in red. Note that for the a given DM mass, the Higgs-DM coupling in case G is much 
larger than in case H. 


Since 6a and 5c are related to the quartic parameters Aj’s, they are constrained from 
unitarity bounds and are required to be small. However, regardless of any limits on Aj from 
unitarity, relic density studies show that both 6a and 5c must be relatively small to allow 
for co-annihilation between particles which is crucial in the heavy DM mass region. The 
upper limit follows the following rough rule: co-annihilation effects take place when the mass 
difference between co-annihilating particles is of the order of 20% of their mass. 

The lower bound on 6a comes from direct detection experiments where a degeneracy 
between Hi and Ai leads to the scattering through the Z boson which is tightly constrained. 
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Further, the above limits in Eq. (23) are in agreement with LEP searches for exotic particles. 
Finally, below 5a,c ~ 0.1 GeV there is no visible difference in the results. 

As mentioned before, a large difference between 5a and 5c violates relic density limits. 
In the cases we study below, we have set 5a = 5c = 5 for simphcitj|^ The mass splitting 
between the two generations. A, is proportional to and therefore unconstrained by unitarity. 
The maximum value for A is proportional to the arbitrary maximum value chosen for pf2- 
In general, we allow for A to vary in the following region 


100 keV < A < 200 GeV. (24) 

For large A values, A > 20 — 50 GeV (depending on ttihi, since the story-changing mass 
splitting is roughly 20%mii/J, co-annihilation effects between the two generations are not 
strong enough to compete with the standard (co)annhilation between Hi, Ai,Hf, in which 
case the I(2-|-1)IIDM acts just like the I(1-|-1)IIDM. The second generation is effectively 
decoupled from the hrst generation and does not influence DM phenomenology. Therefore, 
scenario H is realised for: 

A > 20 GeV scenario H (25) 

The exact value of A, when above ~ 20 — 50 GeV does not make any signihcant difference 
in the relic density calculations. 

For small values of A, the co-annihilation effects between all particles are important. The 
smaller A is, the more relevant particles from the second generation are for DM studies: 
for A 1.5 GeV the relative contribution to relic density calculation coming from particles 
from the lighter generation to the heavier generation is 70%-30%. For A 0.0001 GeV this 
relation is 50%-50%. The case G, is therefore realised when A varies in the following window 


0.0001 GeV < A < 20 GeV scenario G (26) 

The closer A gets to this upper limit, the weaker the coannnihilation effects and the more 
scenario H is realised. Finally, notice that, for our studies, the Higgs-DM coupling, ghHiH^, 
is kept within the \ghHiHi \ < 1 range. 

Fig. i illustrates the effect of changing A on the relic density. In all four plots ttidm 
has been set to 400 GeV as the value of A changes, 0.1 GeV in the top left plot, 1 GeV 
in the top right plot, 5 GeV in the bottom left plot and 10 GeV in the bottom right plot. 
In each plot different colours represent different 5s. For small values of A (0.1 GeV), the 
H 1 -H 2 co-annihilation leads to viable relic density values even for large 5 (i.e. the Hi,Ai, 
co-annihilation is absent). For large values of A (10 GeV) H 1 -H 2 co-annihilation does not 
exist and even small values of 5 cannot compensate this lack, thus, the relic density is below 
the acceptable limit. 

In Fig. the value of 5 is set to 0.5 GeV for mcM = 400 GeV in the left plot and 
fncM = 550 GeV in the right plot. In each plot, different colours correspond to changing As. 

®Cases with 5a A 5c do not lead to any new phenomenology and in fact the region of validity for gniHih 
decreases as we increase the difference between 5a and Sc- 


13 




mH,=400 GeV, A=0.1 GeV 


mH,=400 GeV, A=1 GeV 


-d= 0.015 GeV 





mH,=400 GeV, A=5 GeV 

Odm^^ 



. ^ = 0. 
. ^ = 0. 
. ^ = 0. 
= 1 
6 = 2 
• 6 = 2 
6 = A 
• 6 = S 

-^ = 0 . 
- ^ = 0 . 
1 

■ 6 = 2 
6 = 2 
6 = A 
6 = 5 


015 GeV 
1 GeV 
5 GeV 
GeV 
GeV 
GeV 
GeV 
GeV 

015 GeV 

5 GeV 

GeV 

GeV 

GeV 

GeV 

GeV 


Clnuh'^ 



Clnuh'^ 


inH,=400 GeV, A=10GeV 



- 6 = 0.5 GeV 

- 6= 1 GeV 

- 6 = 2 GeV 

.^ = 3 GeV 

- S = 4 GeV 

. 6=5 GeV 

-^=0.015 GeV 

- 6 = 0.5 GeV 

- 6= 1 GeV 

- 6 = 2 GeV 

-^ = 3 GeV 

- S = 4 GeV 

. 6=5 GeV 


Figure 5: The effect of changing A on the relic density. In all four plots ttidm has been set to 
400 GeV as the value of A changes, 0.1 GeV in the top left plot, 1 GeV in the top right plot, 
5 GeV in the bottom left plot and 10 GeV in the bottom right plot. In each plot different 
colours represent different 6s. For small values of A (0.1 GeV), the H 1 -H 2 co-annihilation 
leads to viable relic density values even for large 6 (i.e. the co-annihilation is 

absent). For large values of A (10 GeV) H 1 -H 2 co-annihilation does not exist and even small 
values of 6 cannot compensate this lack, thus, the relic density is below the acceptable limit. 


Note that in the left plot only small values of A lead to viable relic density values, which is 
where case G is realised. In the right plot, small values of A correspond to case G and large 
values of A correspond to case H, and they all lead to acceptable relic density values. Note 
that for A > 50 GeV all curves correspond to the same value. 


4.2 Changes in tudm 

Here we describe several sub-regimes where the DM mass can vary with important charac¬ 
teristics. 

• In the region 

ruDM < 360 GeV (27) 

neither scenario H nor G results in viable relic density value^ This lower limit can be 
reached in case G by very specific points in the parameter space: (a) when the mass 

®In extensions of the I(2+1)HDM with more inert doublets, this tudm limit could be lowered as more 
co-annihilation channels are present due to the extended number of inert particles. 
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Figure 6: Relic density plot for S = 0.5 GeV for ttidm = 400 GeV in the left plot and 
'moM = 550 GeV in the right plot. In each plot, different colours correspond to changing As. 
Note that in the left plot only small values of A lead to viable relic density values, which is 
where case G is realised. In the right plot, small values of A correspond to case G and large 
values of A correspond to case H, and they all lead to acceptable relic density values. Note 
that for A > 50 GeV all curves correspond to the same value. 


splitting between all particles is tiny and all particles are almost degenerate in mass 
(up to (7(100 keV) mass splitting to avoid direct detection limits), (b) when ghHiH^ is 
close to 0. 

• In the region 

360 GeV < tudm < 525 GeV (28) 

only scenario G leads to acceptable values of DM relic density for specihc values of A 
in the 0.0001 GeV < A < 20 GeV and \gH^Hih\ ^ 0.3 window. As a rule of thumb, 
smaller A allows for a wider viable region in the parameter space and the larger mui 
is the larger \ghH-i_Hi \ must be. 

• In the region 

535 GeV < tudm <1-5-2 TeV (29) 

both scenarios, G and H, lead to viable values of DM relic density. The appropriate 
value of gHiHih coupling depends on the DM mass in each case. In scenario G, couplings 
are generally larger compared to scenario H and a larger DM mass requires larger values 

of the gH,H,h coupling. 

Fig. 0 is meant to summarise all that was said above in one plot with Gi and Hz repre¬ 
senting certain points in the parameter space corresponding to cases G and H, respectively. 
The shaded region is where the I(2-|-1)IIDM has acceptable relic density results. To the left 
of the vertical dashed line case G is realised and to the right of it both cases G and H are 
realised. Generally the outermost parts of the shaded region are populated by Gi since they 
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correspond to larger Higgs-DM couplings whereas the innermost parts of the region corre¬ 
spond to case H. However, depending on the values of S the Higgs-DM coupling changes and 
Gi points can appear close to the QhHiHi = 0 line as well, which is apparent in comparing 
the points G5, HI and G6. Other points are shown on the plot for different values of A and 
6 for comparison. 
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Figure 7: Here Gi and Hi represent certain points in the parameter space corresponding to 
cases G and H, respectively. The shaded region is where the I(2-|-1)HDM has acceptable 
relic density results. To the left of the vertical dashed line case G is realised and to the right 
of it both cases G and H are realised. Generally the outermost parts of the shaded region 
are populated by Gi since they correspond to larger Higgs-DM couplings, and the innermost 
parts of the region correspond to case H. However, depending on the values of 6 the Higgs- 
DM coupling changes and Gi points can appear close to the ghHiHi = 0 line as well, which 
is apparent in comparing the points G5, HI and G6. Other points are shown on the plot for 
different values of A and 6 for comparison. 


4.3 Direct detection 

Direct detection experiments, which are mostly designed to hunt for the standard EW-scale 
WIMP, are the most sensitive to DM masses of the order of 100 GeV. For heavier DM masses 
the sensitivity of these experiments drops signihcantly. The most recent LUX results set the 
limit of the DM-nucleon scattering cross-section to be udm-n ~ 10“® pb for DM masses 
^ 500 - 1000 GeV [I?]. 

In the I(2-|-1)HDM, similar to the other scalar DM models, DM candidate can be detected 
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through elastic scattering on nuclei through the exchange of a Higgs particle. Therefore, 
o'DM-N will depend on the value of DM-Higgs coupling, and the DM mass. 

These results are presented in Fig. |8] where the shaded region corresponds to the probed 
phase space of the I(2+1)HDM for various choices of A and 6 (as shown in Fig. [^, all 
of which have relic density in agreement with Planck measurements. Also, results for the 
benchmark points studied in section are presented explicitly. We found that the current 
experimental limits do not constrain the heavy DM mass region, neither for case H nor for 
case G, even though the Higgs-DM coupling is larger in the latter case. 

Recall that for certain choices of A and 6 one can obtain proper DM relic density for 
OHiHih ~ 0. This leads to a strongly suppressed scattering cross-section, which may not be 
detected as it lies within the coherent neutrino-nucleus scattering regime |18j . 

Fig, I also shows a limit from the future XENONIT experiment [19], with a proposed 
sensitivity of the order of 10“® —10“^° pb (dashed black line). We expect the next generation 
of DM detectors, such as XENONIT, to be able to test a large portion of the parameter 
space of the heavy DM in the I(2-|-1)HDM for < 1 TeV. 

Cf(DM.N} [pb] 

- LUX 

-XENON IT 

- V background 

• case G 

• case H 

mn, [GeV] 

400 500 600 700 800 

Figure 8: The DM-nucleon scattering cross-section for the I(2-|-1)HDM in comparison with 
the direct detection limits set by LUX (black line) and projected XENON IT (black dashed 
line). Coherent neutrino scattering limit is also shown (dashed blue line). The shaded region 
corresponds to points with relic density in agreement with Planck measurements; results for 
benchmark points for case G (A = 1 GeV, 5 = 1 GeV, blue) and case H (A = 100 GeV, S = 
1 GeV) are presented. 

As a hnal note to this subsection, we should like to mention that a viable intermediate 
DM mass region, ^ 160 GeV, regarding relic density studies, has been found 

in the I(1-|-1)HDM. The correct relic density in this region is obtained due to cancellations 
between different diagrams contributing to DM annihilation into gauge bosons {W^W~ and 
ZZ). In [20] it was shown that this scenario is realised if the inert particles, in particular the 
charged scalar, are heavy enough, ~ 300 — 500 GeV. A relatively large DM-Higgs coupling 
is also required for the DM in this mass region to stay within relic density limits, however. 
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this large Higgs-DM coupling is excluded by LUX. Similarly, we did not find any solutions 
in the medium mass region with viable relic density and in agreement with direct detection 
experiments in the I(2+1)HDM. 


5 Heavy Higgs DM at the LHC in the I(2+1)HDM 

A scalar DM candidate is a stable particle with limited interactions with all SM particles and 
therefore it cannot be directly detected at the LHC. However, its presence can influence the 
detectable properties of SM particles. One way to ascertain the influence of DM candidate 
on the properties of a Higgs particle is to look at the Higgs invisible decays, h —)■ SS^ 
where S' is a scalar DM candidate with mass below mh/‘2. Invisible decays of the SM-like 
Higgs particle in the I(2+1)HDM were studied in [T^ 1^ . where limits for the mass of a 
light DM candidate combined with Planck limits for the relic density measurements provided 
constraints comparable or stronger than those from direct detection experiments. 

Another strategy, useful for a heavy DM particle, is to look for a high pt monojet or a two 
jet/two lepton signal, accompanied by a large missing transverse energy The mono jet 
signature in the I(2+1)HDM, pp —)■ HiHi + jet, corresponds to h coupling to an invisible 
pair of DM particles (yielding the large ^t) with produced in association with an energetic 
quark or gluon jet. The following processes are considered in our analysis. 


1. gg —)■ gHiHi (Fig. 16) via a triple gluon and a hgg effective vertex. Note, that the hgg 
effective vertex in the I(2+1)HDM is the same as in the SM, as the Higgs production 
here is not modified by presence of additional scalar states. This is the dominant 
contribution to the monojet process, as the gluon fusion is an enhanced production 
mechanism for the Higgs particle. 

2. qq ^ gHiHi (Fig. ID , where q = u,d, c, s, b. The dominant contribution comes from 
the s-channel via the gqq tree-level vertex and the hgg effective coupling (Fig. 17r). 


3. qg —)■ qHiHi (Fig. 18), where q = u,d,c,s,b. The dominant contributions here come 


from gb —)■ HiHib with the Higgs boson radiated off of the b quark legs (Fig. 18 r) and 


qg —)■ qHiHi t-channel via a gqq tree-level vertex and the hgg effective coupling (Fig. 

Pb). 


Note, that all above processes contain the h — HiHi vertex, therefore a strong dependency 
on the guiHih coupling is expected, i.e., a significant difference between scenarios G and H, 
as discussed in the previous section. 

For the studies of pp —)■ HiHi + 2jets we have considered the Vector Boson Fusion (VBF) 
process of the form HiHig^qi, with q = u,d where a pair of DM particles (with large 

^t) is produced by the neutral (Fig. 19) or charged (Fig. 20) VBF processes, either directly 


or mediated by the Higgs particle or another neutral scalar. 
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We have also considered the Higgs-Strahlung (HS) processes of the form qiqj —)■ V*HiHi 


where a pair of DM particles is radiated off the Z or boson leg (Fig. 21 and Fig. 22 
respectively), either directly or mediated by the Higgs particle or another neutral scalar. 

Notice, that in the dijet searches only one diagram out of each set depends on the gniHih, 
therefore we expect smaller differences between scenarios H and G than in the monojet 
searches. The strength of the other diagrams is set by the gauge interactions. 

Also, it is important to stress that, given our initial choice of parameters, i.e. introducing 


the k = 1 relation between the doublets (see Eq.(12)), we have limited the number of possible 


diagrams, because vertices of the type ZHiAj and W^HiHj {i ^ j) are absent when k = 1. 
Relaxing this initial assumption would in principle not only influence the evolution of DM 
relic density, but could also lead to a possibly stronger difference between scenarios G and H 
in the dijet analysis. 

In the following subsections we present results for the mono jet and dijet analysis, for the 
14 TeV LHG. The following selections were used. 

1. For the monojet searches, we require the following cuts on the transverse momentum 
of the jet, and the pseudo-rapidity of the jet, pG 


plp > 20 GeV and \r]^\ > 2.5 


(30) 


2. For the dijet searches, we require the following cuts on the invariant mass of the two jets, 
and the difference between the pseudo-rapidity of the forward and backward 
jet, 

M{j,j) > 700 GeV and Ivj — vll >4 (31) 

Galculations were done with the aid of LanHEP [22] and GalcHEP [23] packages. 


5.1 Mono jet results 

In Fig. [^results for monojet signals of scenarios G (5 a = 5c = 1 GeV, A = 1 GeV) and H 
(5a = 5c = 1 GeV, A = 100 GeV) are shown. For comparison, we also present results for the 
I(1-I-1)HDM, with 5a = 5c = 1 from [12]. The DM-Higgs coupling (dehned as A345 in [T2] j 
is the same as gniHih in scenario H for equal DM masses, therefore the monojet diagrams in 
case H and in I(1-|-1)HDM are identical. 

Scenario G, which corresponds to much larger Higgs-DM couplings compared to that of 
scenario H or the I(1-|-1)HDM, results in a signihcantly larger cross-section in the monojet 
process. Also the special features of the model are more visible in this process. For masses up 
to niHi ~ 450 GeV we observe a rise in the cross-section connected to an opening of the phase 
space combined with an increasing Higgs-DM coupling. After that peak, the cross-section 
decreases with increasing DM mass regardless of the rising of gtHiHi- 

Notice that for the lower masses, the difference between scenario G and H is signihcant, 
as every diagram involved in the monojet process contains the hHiHi vertex, whose coupling 
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differs significantly in cases G and H. Notice the low end of the allowed mass region in case 
G with a large cross-section for the mass region which is not even accessible by case H or the 
I(1-|-1)HDM. As the DM mass grows, results for both cases get closer together, stabilising 
for the very heavy mass region in the decoupling limit. 


a[pb] 


Monojet 



mH,[GeV] 


— case G 

— case H 

• the IDM 


Figure 9: Mono jet searches for cases G and H. For comparison we also provide mono jet 
searches for the I(1-|-1)HDM which resemble case H in the I(2-|-1)HDM as expected. 


5.2 Dijet VBF results 


Fig. [^presents values of the dijet cross-section for scenarios H and G in terms of the DM 
mass. The difference between cases H and G is less prominent compared to the monojet 
analysis, as only one of the involved diagrams in this process depends on the value of the 
Higgs-DM coupling (see Figs. IT^ and 20 for diagrams involved in the neutral and charged 


VBF processes, respectively). We can still observe some differences in the lower mass range. 
The cross-section for case G is generally larger (as is the QhHiHi strength) than in case H. 
Also, charged channels have slightly larger cross-sections than the neutral ones, since the 
cross-section for producing the boson is larger than the cross-section for producing the 
Z boson. For the heavier masses all results, for both scenarios G and H, as well as for charged 
and neutral channels, tend to converge. 


5.3 HS results 

HS signatures depend on the W and Z decay patterns. While at the LHG, leptonic signatures 
are preferred, hadronic ones are also possible. The latter potentially interfere with the VBF 
topologies, but the effect is small so that we can safely ignore it here. 

The results of the (on-shell) HS process cross-sections in terms of the DM mass are 
presented in Fig. It is clear that the general picture is different from the VBF studies. 
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Figure 10: The cross-section in dijet VBF processes in terms of the DM mass in charged 
(top) and neutral (bottom) processes. 


The largest cross-section again appears in the lower mass region where only case G can be 
realised. Similarly to the VBF case, the charged channels have larger cross-section compared 
to the neutral channels since the production in theses processes has a larger cross-section 
than the Z boson production. All cross-sections in neutral and charged processes in both 
cases G and H converge to a similar value as tudm approaches very large values. 

Note that the cross-section in both G and H are similar in the region of DM mass above 
525 GeV, i.e., where both cases can be realised (unlike in the VBF processes). This similarity 
is the result of the fact that the difference in the Higgs-DM coupling does not translate into a 
difference in the cross-section between cases G and H. To explain this similarity let us focus, 
e.g., on the neutral VBF and HS processes (Fig. 19 and Fig. imQ 


In the neutral VBF process, out of all the involved diagrams (Fig. 19i,b,c) there is only 


one diagram. Fig. 191, which depends on the Higgs-DM coupling. The cross-section of this 


^The same argument applies to comparing the charged VBF and HS processes. 
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Figure 11: The cross-section in HS processes in terms of the DM mass in charged (top) and 
neutral (bottom) processes. 


diagram (ah) for a given ttidm relative to the cross-section of all three diagrams involved 
(cTtot) is (Jh/cTtot = 0.1445 for case G and Uh/utot = 0.2416 for case H. Recall that the main 
difference between cases G and H is that for a chosen muM the Higgs-DM is larger in case G 
than in case H. We conclude that the Higgs-mediated diagram plays a much more important 
role in case H than it does in case G. As a result, even though the gtHiHi coupling is much 
smaller in case H, the total cross-section does not fall far below the total cross-section in case 
G, which is depicted in Fig. [Toj 


Now, let us consider the HS neutral processes (Fig. 21i,b,c). Again, only one diagram. 


Fig. 211, depends on the Higgs-DM coupling. Repeating the procedure above, we calculate 


the relative cross-section of this one diagram (ah) relative to the cross-section of all diagrams 
involved {(Jtot)- For a given moM, we obtain <Jhl<Jtot = 0.9736 in case G and ahjotot = 0.9708 
in case H. So, this diagram plays only a slight role in case H compared to case G. We 
therefore conclude that the difference in ghHiHi coupling does not play an important role in 
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distinguishing case G and H in the HS processes. Thus, we do not expect to see a difference 
between cases G and H, which is depicted in Fig. 11, with the important exception that 


case G allows for the heavy Higgs DM mass to be below 525 GeV, making its cross-section 
significantly larger. 


6 Conclusion 

We have calculated the relic density for heavy Higgs DM in the I(2-|-1)HDM and shown 
that the prospects for its discovery at both DD experiments and the LHG are significantly 
enhanced as compared to the I(1-|-1)HDM, where the heavy Higgs DM particle must have 
a mass above about 525 GeV and is weakly coupled to the observed Higgs boson. Adding 
a second inert Higgs doublet helps to make the heavy Higgs DM region accessible to both 
DD and the LHG, by either increasing its couplings to the observed Higgs or lowering its 
mass to 360 GeV < moM, or both. In particular we have presented LHG signatures of the 
I(2-|-1)HDM in the monojet, VBF (dijet) and HS processes and shown that the prospects for 
heavy Higgs DM discovery are significantly brighter for all channels. 

In DD experiments, although the standard values of annihilation cross-section for the 
heavy Higgs DM masses in the I(2-|-1)HDM are well below current experimental exclusion 
limits for DM decaying into pairs of gauge bosons or fermions, the prospects for a future DD 
discovery remain open due to the complementary nature of collider vs cosmological limits 
and the fact that the DD cross-sections are higher than in the I(1-|-1)HDM. 

Turning to indirect detection signatures of the I(2-|-1)HDM, there is the possibility of in¬ 
ternal bremsstrahlung in the processes of HiHi —)■ generated through the exchange 

of any of the two charged scalars shown that one can expect such signatures in 

the I(1-|-1)HDM [21], mediated by a charged scalar in the t-channel, which would correspond 
to scenario H considered in the I(2-|-1)HDM. In principle, the signal could even be stronger 
for scenario G, as the scalar couplings are enhanced. 

Finally we comment on the observed h —?• 77 channel where, in the I(1-|-1)HDM, only 
in the heavy DM mass region are both proper DM relic density and (minimal) enhancement 
in the h —?• 77 channel realised. By contrast, in the I(2-|-1)HDM there exist two charged 
scalars, Hf and which contribute to the h —)■ 77 loop which may enhance the rate for 
a wide range of parameters. 

In conclusion, adding a second inert Higgs doublet significantly improves the prospects 
for observability of heavy Higgs dark matter in future experiments both underground and at 
the GERN LHG. 
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A Feynman rules in the simplified I(2+1)HDM 

Here we present the Feynman rules of the model. 

H+H^h, H+H^h X23V 

HiHih, H 2 H 2 h (A23 + A23 + 2A2)- 

AiAih, ^2^2^ (-^23 + -^23 ~ 2A2)- 

HtH^l Anew + g' cos ew){K + K'Y 

H+H^Z, H+H{Z ^ {g cos Ow - g' sin Ow) {K + K'Y 

HiH2W^ 'j-cosYh-eYiK + K'Y 

HfH2W^ y sin(0;, - + K'Y 

HfAiW^, H^A2W^ I cosiOa - OYiK + K'Y 

H^A,W^, HfA2W^ I sin(0, - dY{K + K'Y 

HiAi, H 2 A 2 Z ^{g cos Ow + g'An ew)cos{eh-Oa){K + K'Y 

H 2 A 1 , H 1 A 2 Z ^{g cos 9w + g'sin 9w) An{9h - 9a){K + K'Y 

where K and K' are the momenta of the associated particles in the decay channel and dw 
is the Weak mixing angle. The Yukawa couplings in the model are identical to those of the 
SM. 

B Feynman diagrams for relic density calculations 

Here we present the DM (co)annihilation diagrams which play a role in our relic density 
studies. 
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Figure 12: For light DM masses the most important channel for the annihilation of DM 
particles is the Higgs-mediated process (a). Coannhilation with other neutral scalars could 
have a signihcant effect on the relic density (b). 




Figure 13: For heavy DM masses, the processes involving one or two virtual gauge bosons, 
also contribute to the total annihilation cross-section and affect the relic density. 



Figure 14: 


Heavy DM (co)annihilation diagrams with pure gauge boson hnal states. 
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Figure 15: Heavy DM (co)annhilation channels involving the SM-like Higgs boson. 


C Feynman diagrams for the LHC analysis 

Here we present the DM (co)annihilation diagrams which play a role in our LHC studies. 

C.l Diagrams with monojet final states 



Figure 16: Relevant monojet diagrams with initial gluon states {gg hg ^ gHiHi) 
containing triple gluon vertex and an effective ggh vertex. 
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Figure 17: Relevant monojet diagrams with initial quark states {qq —)■ gHiHi + diagrams 
with initial particles reversed) containing ggh effective vertex, where q = u,d, c, s, b. 



(a) (b) (c) 

Figure 18: Relevant monojet diagrams with initial quark and gluon states {qg —)■ qHiHi + 
equivalent qg —)■ qHiHi diagrams + diagrams with initial particles reversed) containing ggh 
effective vertex, where q = u,d, c, s, b. 
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C.2 VBF diagrams with dijet final states 







Qi 





(b) 


Figure 19: Relevant VBF diagrams with dijet final states {qiqj —)■ HiHiqiqj) 
with only neutral intermediate gauge bosons, where q = u,d. Note that 
only one of the involved diagrams in this process depends on the value of 
the Higgs-DM coupling. 



(b) 


(c) 


Figure 20: Relevant VBF diagrams with dijet final states —)■ HiRiq^qi) 

with only charged intermediate gauge bosons, where q = u,d. Note that 
only one of the involved diagrams in this process depends on the value of 
the Higgs-DM coupling. 









C.3 HS diagrams with (on-shell) gauge boson final states 



(a) (b) (c) 

Figure 21: Relevant HS diagrams with (on-shell) neutral gauge boson hnal 
states {qiqi —)■ HiHiZ*) and only neutral intermediate gauge bosons, where 
q = u,d. Note that only one of the involved diagrams in this process 
depends on the value of the Higgs-DM coupling. 



Figure 22: Relevant HS diagrams with (on-shell) charged gauge boson hnal 
states {qiqj —)■ HiHiW*^) and only charged intermediate gauge bosons, 
where q = u,d. Note that only one of the involved diagrams in this process 
depends on the value of the Higgs-DM coupling. 
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